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1 Introduction 



In spite of its long history, statistical mechanics for nonequilibrium systems remains to 
be established. The central principle for the general theory is still lacking. Recently, a 
lot of analytical and numerical studies of specific nonequilibrium models have revealed 
their highly interesting behaviors. Among them, exactly solvable models are of great value 
because physical quantities and their relations are obtained explicitly PQEIISI- 

The one- dimensional asymmetric simple exclusion process (ASEP) is an exactly solvable 
stochastic process of many particles jSHUElin]- Particles subject to the exclusion interaction 
perform random walk with different rates in the right and the left direction. Even after the 
relaxation to the steady state, the ASEP has nonvanishing particle current, maintained 
in a nonequilibrium situation. The ASEP has become one of the standard models for the 
study of nonequilibrium statistical mechanics because of its simplicity, rich behaviors and 
wide applicability. For instance, the ASEP is regarded as a primitive model of kinetics 
of biopolymerizationjTj, traffic flowjHI, formation of shocks[9j. It is a discrete version of 
a hydrodynamic system obeying noisy Burgers equation [lOj. It also appears in a kind of 
sequence alignment problems in computational biology ^1]. 

In this article we study the ASEP with open boundary conditions. The model has 
in total five parameters 7,5 and q. The parameter q is related to the asymmetry 
of the particle hopping; a (resp. 5) and 7 (resp. /3) represents the particle input (resp. 
output) rate at the left and the right boundary respectively. The stationary properties 
of the model has attracted much attention. Among several reasons, one of them may be 
that it shows a boundary induced phase transition |12j: even the bulk properties of the 
system change drastically depending on the values of the boundary parameters. For the 
case where particles hop only in one direction (g = 7 = 5 = 0), the current and the density 
in the thermodynamic limit was calculated in ^3 E]- In QB]; the authors introduced 
a novel algebraic formulation of the problem, which we call the matrix method in the 
following. This is a method of constructing the exact stationary measure of the model in 
terms of matrix products. It has become a standard techniques for investigating stationary 
properties of the ASEP and related models. A lot of generalizations have been considered 
including discrete time models and multispecies models. 

In ^^-i for the continuous time setting, the algebraic formulation of the solution has 
been given for case with the full five parameters although the computation of physical 
quantities in the thermodynamic limit was performed only for the case where q = J = 
5 = 0. In jTH] the five parameter case was treated. The current in the thermodynamic 
limit was calculated and the phase diagram was obtained. The analysis, however, included 
some approximations. In addition, the computations of other quantities such as the density 
using the same direction of argument seemed difficult. In [IHIIIZI, the case where 7 = 5 = 
was studied. The analysis of was based on the findings that the stationary state of the 
ASEP for the case is related to the g-orthogonal polynomials called the Al-Salam-Chihara 
polynomials. This method has also admitted the computation of the density profile in JH] • 

The objective of this paper is to investigate the most general case with the five pa- 
rameters and further clarify the intimate relationship between the stationary state of the 
ASEP and the theory of g-orthogonal polynomials. We exactly calculate the partition 
function, the n-point functions and other physical quantities. Each quantity is written in 
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an explicit integral form in terms of the weight function of the Askey- Wilson polynomials. 
The asymptotic behaviors in the thermodynamic limit are easily evaluated using analytical 
methods for integrals. As a result, we confirm the phase diagram for the current obtained 
in [H!. 

The paper is organized as follows. In section |^ the definition of the model and physical 
quantities of the ASEP are introduced. The matrix method is also explained. In section 
OJ properties of the Askey- Wilson polynomial are reviewed. The recurrence relation and 
the orthogonality relation play key roles in deriving the integral formula of the partition 
function. In section |3J two representations of the algebraic relations of matrices are given. 
One of them is related to the g-Hermite polynomial, which is the special case of the Askey- 
Wilson polynomial. It is simpler but the parameters are restricted. The other is related 
to the Askey- Wilson polynomial with full five parameters where no restriction is needed. 
The connection between the two representations is also mentioned. Section El deals with 
the case where q = 1- In section [HI integral formulae of the partition function are derived. 
In the thermodynamic limit, the asymptotic behaviors of the partition function and bulk 
quantities of the ASEP are examined. There appears nonequilibrium phase transition. In 
section [71 using formulae in the g-calculus, integral formulae of the n-point functions are 
derived. The last section is devoted to the conclusion. 

2 Model and Matrix Method 

The one-dimensional ASEP is a model of randomly hopping particles with exclusion in- 
teraction on a one-dimensional lattice [H 121 El lH El ISj- Let us denote the system size 
as L. Each site z = 1, ■ ■ ■ , L of the lattice can afford at most only one particle; two or 
more particles cannot stay on the same site due to the hard-core exclusion interaction. A 
particle hops to the nearest right (resp. left) site with the probability pndt (resp. pidt) 
during an infinitesimal time dt if the target site is empty. By rescaling time, we can set 
Pr = 1 and Pl = Q- In the case of open boundaries, a particle enters the system at the 
left (resp. right) boundary with rate a (resp. S) if the site 1 (resp. L) is empty, and is 
removed at the left (resp. right) boundary with rate 7 (resp. j3) if the site 1 (resp. L) 
is occupied. Figure ^ illustrates how the system evolves. Obviously, we can think of the 
reflection symmetry in this finite system: 

site — 2 + 1, i = 1, ■ ■ ■ , L 

rate q ^ q^^, a ^ q '^5, {3 ^ q^'^'J, 7 ^ q^^P, S ^ q^^a. (2.1) 

Therefore, the case where g > 1 is transformed into the case where g < 1 by the symmetry. 
In the following, we only assume the boundary parameters a,/3, 7,5 are strictly positive 
unless otherwise stated. 

The configuration of the system is labeled by {ti,T2,--- ,tl) where = 0, 1 is the 
particle number at site i. Let P{ti,T2, ■ ■ ■ ,TL]t) denote the probability in the configuration 
(ti,T2,-- - ,tl) at time t. The time evolution of P(ri,r2,--- ,TL;t) is described by the 
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master equation, 

d , 

— P{Ti,T2, ■■■ ,TL]t) = 

at 

5^,,i[aP(0, rs, ■ ■ ■ , tl] t) - 7P(1, ra, ■ ■ ■ , r^; t)] 
+^ri,o[7^(l, ^2, ■ ■ ■ , 'Tl\ t) - aP(0, r2, ■ ■ ■ , r^; t)] 

L-l , 

+ ^(r„r.+i)=(i,o)[g^'(- ■ ■ ,ri = 0,ri+i = 1, ■ ■ ■ ;t) - P(- ■ ■ ,ri = l,ri+i = 0, ■ ■ ■ ; t)] 
i=i 

+(5(^„r,+i)=(o,i)[^(--- >T-i = l,ri+i = 0, ■■■ ;t) -qP{--- ,Ti = 0,Ti+i = l,--- 
+6r,,o[f3P{T,, ■■■ , Tl_,, 1; t) - 6P{n, ■■■ , Tl_,, 0; t)] 

+5.,,i[5P(ri, ■ ■ ■ , TL-i, 0; t) - (3P{n, ■ ■ ■ , r^.i, 1; t)]. (2.2) 

Hereafter, we focus on the stationary situation of the ASEP. According to the matrix 
method |13j, the stationary solution P(ri, r2, ■ ■ ■ , tl) for the master equation is obtained 
in the form, 

L 

P{ru r2, ■ ■ ■ , tl) = ^{W\ [^(r.D + (1 - r,)E)|V). (2.3) 



i=l 



Here D and E are matrices; D (resp. E) corresponds to an occupation (resp. emptiness) of 
a particle. {W\ and \V) are vectors; {W\ (resp. \V)) corresponds to the left (resp. right) 
boundary. The normalization constant Zi is given by 

Zi = (iy|C^|\/), (2.4) 

where 

C = D + E. (2.5) 

We call Zi the partition function. One can show that (j2.3p gives the exact steady state of 
the system if the matrices and the vectors satisfy |13j 



DE-gED = D + E, (2.6) 

{W\{a£-^D) = {W\, (2.7) 

{(3D-5£)\V) = \V). (2.8) 

Physical quantities are written in the form of matrix products. The average particle 

number at site z, (xj), where the bracket means the average over the stationary probability 
distribution (j2.3j) . is written as 

(r.) = ^{W\C-'DC'^-^\V), (2.9) 

and the two-point function (rjTj) as 

(r,r,) = ^{W\C~^DC^-'-^DC^-^\V). (2.10) 
Zl 



The n-point functions are expressed similarly. The particle current through the bond 
between the neighboring sites from left to right, which is defined by J = (rj(l — Tj+i) — 
g(l - ri)ri+i), is simply given by 



J 



(2.11^ 



The site-dependence of J is unnecessary since the current is uniform for the steady state. 
One can introduce the fugacity in Zl as 



ZL{e) = {w\{e^+^f\v). 



(2.12) 



The square of fugacity in the definition is just for convenience. Then, the mean and the 
variance of the particle density in the whole system are calculated from the first and the 
second derivative of Zii^"^) 



{p) = Tei^^ogz,ie) 



d 



log^L(e') 



(2.13) 
(2.14) 



As is familiar in equilibrium statistical physics, the bulk quantities are written in terms of 
the partition function. 



3 Askey- Wilson Polynomial 



The Askey- Wilson polynomial is a g-orthogonal polynomial with four free parameters be- 
sides q. It resides on the top of the hierarchy of the one-variable g-orthogonal polynomial 
family in the Askey scheme [201 IS 122] ■ We list several properties of the Askey- Wilson poly- 
nomial needed in the following sections. For the purpose we introduce notations heavily 
used in the g-calculus. In this section, |g| < 1 is assumed. The g-shifted factorial: 



[ai, (32, ■ ■ ■ 5 (^s', 



g)n=nn(i 



(3.1) 



r=l k=0 



The basic hypergeometric function: 



ai. 



k=0 



(ai, ■ ■ ■ , ttr] q)k 
(6i, - ■■ ,bs,q;q)k 



(3.2) 



Then, the Askey- Wilson polynomial Pn{x) = Pn{x; a, b, c, d\q) is explicitly defined by 

^g-", q^'-^abcd, ae*^ ae"^^ 



Pn{x) = a ^{ab, ac, ad; q)n 403 



(3.3) 



a6, ac, ad 

with X = cos6 for n G Z+ := {0, 1, 2, ■ ■ ■ }. It satisfies the three-term recurrence relation, 

AnPn+lix) + BnPn{x) + C„P„_i(x) = 2xP„(x), (3.4) 
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with Po{x) = 1 and P_i(x) = 0, where 
1 - q^-^abcd 

" (1 - q^^-^abcd){l - q^^ahcd) ' ^ ' 

(3.6) 

(1 - g")(l - g"-ia6)(l - g"-iac)(l - g"-iaci)(l - g"~i6c)(l - - g^-^crf) 

(1 - g2n-2a5cd)(l - q^^-^abcd) ' 

(3.7) 



Cn 



and 

s = a + 6 + c + c?, s' = a'^ + b'^ + c'^ + d~^. (3.8) 

From this relation, we see that Pn{x;a,b,c,d\q) is symmetric with respect to a, b, c and 
d. The orthogonahty relation depends on the parameter region. If \a\, \b\, |c|, \d\ < 1, the 
orthogonality reads 

—w{cOs6)Pm{cOsd)Pn{cOsd) = hnSmn, (3.9) 

Ztt 



where 



(e2^^e-2^^;g)oo 



u;(cos6') ^^^.^^ ^^^^^ ^.^^ ^^.g^ , (3.10) 



(3.11) 



/i„ (1 — g" abed) (g, a6, ac, ad, be, bd, ed; q)n 
Hq (1 — q^^~^ abed) [abed; q)n 

^ (abed; q)^ 

(g, ab, ae, ad, be, bd, ed; q)ao 



If we rewrite (13. 9j) with z = e*^, we have 



f i±fl\ p„ f i±£l) p. { i±£^ ) = (3.13) 



where the integral contour C is a closed path which encloses the poles at z = aq'', bq^, eq^, 
dq^ {k e Z,+) and excludes the poles at z = (aq'^)'^, {bq^)~^, {eq^)~^, {dq^)~^ [k G Z+). In 
the other parameter region, the orthogonality is continued analytically. Putting m = n = 0, 
we have the celebrated Askey- Wilson integral, 

dz {z'^,z-^;q)^ _ {abed;q)^ 



Aniz (az, a/z, bz, b/z, ez, ej z, dz, d/ z]q)oo (<?, ob, ae, ad, be, bd, ed; q)oo 

For the case ofa = b = e = d = 0, the Askey- Wilson polynomial is reduced to the 
g-Hermite polynomial, 

Hnix\q) = Pn{x;0,0,0,0\q). (3.15) 



Its recurrence and orthogonality relation are those with a = b = c = d = in the above for- 
mulae for P„(x; a, b,c,d\q). In what follows, we sometimes write briefly the corresponding 
weight function or coefficients as f^'^^ := f\a=b=c=d=o- We introduce the g-binomial, 

^n(^.^) = E . ^"-v, (3-16) 

[Q, (i)k{q, qjn-k 



k=0 



which satisfies the following recurrence relation, 

Fn+iix,y) - ix + y)Fn{x,y) + (1 - = 0. (3.17) 

with Fo{x,y) = 1 and F-i{x,y) = 0. The g-Hermite polynomial is the special case of it: 

H4cose\q) = F4e'',e-''). (3.18) 

As to Fn{x, y), the g-Mehler formula is important: 

Fn{v,w)Fn{x,y) ^ {vwxy^^; q)^^ 
^ (?;g)n {vxi,vyi,wxi,wyi]q)J 

for \vy^\, \wxE,\i \wyS,\ < 1. Combining ()3.17|) with this, we obtain another kind of 

bilinear function: 

y> Fn{v,w)Fn+i{x,y) _ x + y- {v + w)xy^ {vwxy^^;q)oo ^g) 
^-^ {a:a)r, 1 — vwxv^'^ (vxf,vv£,wxf,wv£:a)^^ 



1=0 



for \vx^\, \vy$,\, \wx^\, \wy^\ < 1. 

When abed ^ 0, the symmetry relation for 4^3, 



493 



d,ej 



implies the symmetry relation for the Askey- Wilson polynomials, 

Pn{x]a,b,c,d\q) = {-abcdYq^'''^''-^^'^Pn{x-a-\b-\c-\d-^\q-^). (3.22) 

From this symmetry, it is enough to study the case where |g| < 1. But when abed = 0, this 
relation no longer holds and one has to consider the case where g > 1 separately [221 121] • 
For instance suppose that only one of a, b, c, d is zero. Without loss of generality, we can 
set c? = because of the symmetry in a, b, c, d. The recurrence relation ()3.4|) remains the 
same (with d set to zero), but the orthogonality relation ()3.9p should be modified. P„(a;)'s 
are now orthogonal on the imaginary axis; 

/oo 
du w {sinh u) Prn{i sinh. u) Pn{i sinh. u) = hnSmn, (3.23) 
-00 

where 

, , 1 (iaq'^e^, ~iaq~^e'^,ibq~^e'^, —ibq~^e~^,icq~^ e"^, —icq'^e'"^] q~^) 00 

w{smhu) = ^^-^ , 

logg (-g ^e^",-g g ^)oo 

(3.24) 

hn = (g, ab, ac, be; q)n{q~^, abq~^, acq~^, bcq~^; q~^)oo- (3.25) 

The formulae when more than one of a, b, c, d are zero are similar. For details about the 
case where g > 1, we refer the reader to [IHl El]. The list of formulae exploited in the 
following sections is completed. 
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4 Representation: the case where 



In the subsequent sections we give two infinite-dimensional representations of matrices and 
vectors. One is related to tlie g-Hermite polynomials and the other to the Askey- Wilson 
polynomials. The representations are useful since the g-orthogonal polynomial forms a 
complete set of eigenvectors of the matrix C, which enters in the matrix expression of 
fl2.4|l . This leads to an integral formula of it. In the following the parameters are fixed as 



where 



(4.1) 



o,7 



2a 
1 

2P 



(1 - g - a + 7) ± - g - a + 7)2 + 4^7 
'l-q-(3 + 5)± y^il-q- /3 + 6)^ + A/36 



(4.2) 
(4.3) 



We also fix g as < g < 1 ; the g > 1 case is simply considered from the reflection symmetry 
()2.1|) . or equivalently, 



a,b,c,d,q ^ b \a ^. 



(4.4) 



Since 7,5 are strictly positive, the new parameters satisfy a > 0, 6 > 0, — 1 < c < 
and —1 < (i<OifO< g< 1. The bra and ket notation are used for row and column 
vectors, respectively. In place of matrices D and E, we consider d and e. 



1-g 
C = D + E 



;i+d), 
1 



1 



l-g 



1-g 
(21 + d + e). 



:i + e) 



For them the relations (j2.6|) . (|2.7p and (j2.8j) are replaced by 

de — ged = (1 — g)l, 
{W\[e-{a + c)l + acd] = 0, 
[d - {b + d)l + bde\\V) =0. 



(4.5) 
(4.6) 



(4.7) 
(4.8) 
(4.9) 



4.1 Representation related to the g-Hermite polynomials 

First, we give a representation related to the g-Hermite polynomials. One can check the 
following matrices and vectors satisfy the relation ()4.7p . fl4.8p and ()4.9|) . 



y/T^ 

v^l-g2 























g2 



,(4.10) 
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{W\ = Ua,c [Fo{a, c)/ ^J{q] q)o, Fi{a, c)/ 



\V) = v,,d [Fo{b, d)/^/{q; q)o, F,{b, d)/^{q- g)i, ■ ■ -j , (4.11) 

where Va,c = {a-c; q)^{q', <l)oo^'^- We regard ()4.7|) as the g-deformed commutation relation of 
boson operators [2nil2Sll2Zl 5 and (j4.1Up and (j4.11|) as the g-deformed Fock representation [T7] . 
Note that matrices d and e depend only on q, and vectors {W\ and \V) depend on the 
boundary parameters as well. In this representation, a vector defined by 

{ho{x),hi{x),h2ix),- ■ ■)'^, 



where 



IHx)) 



hn{x) 



(g;g)c 
(?;?)7 



-if„(x|g), 



is eigenvector of the matrix d + e with the eigenvalue 2x: 

(d + e)|/i(x)) = 2x\h{x)). 



(4.12) 



(4.13) 



(4.14) 



Each element of this equation corresponds to the recurrence relation of the g-Hermite 
polynomials. The transposed equation also holds for {h{x)\ := \h{x))^ because d + e is a 
symmetric matrix. The orthogonality relation is rewritten as 



1 



27r 



w^"^(cos^) |/i(cos^))(/i(cose)|. 



(4.15) 



Therefore, the eigenvector of d + e, or at the same time the eigenvector of C from ()4.6p . 
{\h{cos0))]6 G [0,7r]} forms a complete set of vectors in this representation space. Fur- 
thermore, from ()3.20j) with ()3.18j) we obtain the following form factor, 

(g, ^^g; 9)00 (2 cos 0-2^ cos 0) 



(/i(cos0)|A(Od|/i(cos0)) 



(4.16) 



Here a regularization matrix A(^) = diag(l, ^, ■ ' ' ) is introduced. 

In fact, this representation is valid only when |a|, |c|, \d\ < 1, otherwise, the boundary 
vectors are ill-defined. For example, if a > 1, (yV\h{cos9)) becomes a divergent series. See 
the convergence condition for the g-Mehler formula (|3.19|) . In the next subsection, we 
present another representation without such restriction. 



4.2 Representation related to the Askey-Wilson polynomial 

We consider a representation related to the Askey-Wilson polynomials with full five param- 
eters. Here a, 6, c, d are not restricted. Complicated apparently, the followings certainly 
solve the equations (lO) . (Oil and 



dl dl 

di-^ d-^ 

d\ d\ 



-0 ^0 







(4.17) 
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{W\ = 0, 0, ■ ■ ■ ), 1^) = h'J'il, 0, 0, ■ ■ ■ (4.18) 

where 



(1 - q^'^-'^ahcd){l - q^'^ahcd) 

X [6(i(a + c) + (6 + (i)g - a6cc/(6 + rf)g""^ - {6c/(a + c) + a6cc/(6 + d)}q^ 
-bd{a + c)q''+^ + ab^cd\a + c)q^''-^ + ahcd{h + d)q^% (4.19) 



and 



(1 - q^''^'^ahcd){l - q^'^ahcd) 

X [ac(6 + (i) + (a + c)g — abcd{a + c)q^~^ — {ac{b + d) + abcd(a + c)}g" 
-ac{b + d)q''+^ + a%^d{b + rf)g2"-i + abcd{a + c)?^"] , (4.20) 



1 - g"ac " 1 - g"ac ^ ^ 

> - (tbd_j. b ^ A u 221 



;i - q''-^abcd){l - g"+i)(l - g"a6)(l - g"ac)(l - g"arf)(l - g"6c)(l - g"6rf)(l - q^cd)' 
(1 - g2n-ia5crf)(l - q^'^abcdfil - q^^'+^abcd) 

(4.23) 



As before, a vector defined by 

W)) = (po(a;),Pi(x),p2(x),---)^, (4.24) 

where 

Pn{x) = )J^Pn{x] a, b, c, d\q), (4.25) 

becomes the eigenvector of the matrix d + e with the eigenvalue 2x: 

(d + e)\p{x)) = 2x\p{x)) . (4.26) 

Each element of this equation corresponds to the recurrence relation of the Askey- Wilson 
polynomials ()3.4|) . The transposed equation also holds for {p{x)\ := |p(x))^ because d + e 
is a symmetric matrix. The orthogonality relation ()3.9|) is rewritten as 

1 = £ {{z + z-')/2) \p {{z + z-')/2)){p {{z + z-')/2) I, (4.27) 

where C is as described below ()3.13|) . The reflection ()4.4|) corresponds to just taking 
transpositions in the representation. If we set c = = 0, or equivalently 7 = 5 = 0, we find 
the representation is reduced to the one related to the Al-Salam-Chihara polynomials |16j . 
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The relation between the representation here and the previous one is almost obvious. 
They are equivalent if |a|, |c|, \d\ < 1. The connection matrix comes from the transfor- 
mation, 



k 



-1/2 



-\p{cos9)) = G[l'^f^ \h{cos9)). 



We can calculate the elements as 

(ab, ac, ad; q)n 



(4.28) 



n,k a"(ac; q)oo^/K{q] q)k{q] q)oo ^ {q,ab,ad;q) 
The inverse matrix is 



E '''A'":'°"\'''''' fl('''f.'^). (4.29) 



^(a,c) 



(a,c) 



(4.30) 



To prove this, compare the two orthogonality relations ()4.15p and ()4.27p . Replacing 
(a, c) ^ (6, d) and transposing in ()4.28|) . we have 



h, 



-1/2 



■(j9(cos0)| = {h{cos9)\ 



(4.31) 



Then, 



1 = / — w{cos9) |p(cos6'))(]9(cos( 

io 27r 











^w(°)(cos0) |/i(cos0))(/i(cos( 



(a,c) 



1 T 



^(a,c) 
^{b,d) 



I Hermite 



This matrix actually connects the two representations since, for example, {W\ 

Looking at the form factor in the g-Hermite representation ()4.16|) . a nontrivial formula 
results in the Askey- Wilson representation; 



(p(cos0)|e(Odb(cos0)) = 

(2 cos - 2^ cos 9) (g, (^q, ae''^, ae"^"^, ce*-^, ce-*-^, 6e*^ be"^ , de'\ rfe"^^; q). 



where 



0(0 = G'|:>(0 



G 



(a,c) 
{b,d) 



(4.32) 



(4.33) 



The matrix elements are explicitly calculated as 



[e(0] 



1 {abcdC,"^; q)oo{.oh, ac, ad; q)n{ab, be, bd; q)^ 



a-^'b 



nu-k 



"'^ Vhnhk {q, oh^, ac, ad^, bc^, bd, cd^; g)oo 

q'^, q^~^abcd, ab^q^ , bcC, 



•sr^ (q ^abcd,ab^,ad^;q)jq^ 
(g, ab, ad, abcd^'^; q)j 



4(^3 



ab, be, abcd^q^ 



(4.34) 
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A remark on a special case is in order. When we specialize to the case where ab = q^~^ 
{N = 1, 2, ■ ■ ■ ), the infinite-dimensional representation can be reduced to an A^-dimensional 
one. The same is true for other pairwise products of a, b, c, d because of the symmetry in 
a, 6, c, d. The matrices can be taken from the first N x N entry of those of the infinite- 
dimensional matrices, and similarly the first N entry for the vectors: 



d^ d^ 

■•• 







"iV-2 

d^ d^ 







, (4.35) 



(l^(^)| = /.f (l,0,---,0), |yW) = /^f (l,0,---,Or. 



(4.36) 



In this case, the corresponding Askey- Wilson polynomials are called the g-Racah polyno- 
mials. The orthogonality relation is rewritten in N residue sum from the integral form. 
This special case may be useful for trial calculations of many quantities of the ASEP 
since a generally hard task of integration is replaced by finite summation. In fact, this 
finite-dimensional representation is equivalent to the one in [2H1 121] ■ The similarity trans- 
formation is realized by an x matrix U : 



_ {ab,ac,ad;q)n{ab,bc;q)k 



for = 0, 



AT 



- 1. Then, 



q abed, abq 

ab, ad 



(4.37) 



6 
bq ■-. 

6g2 



6-1 
1 






(bq)-' 
1 



bq 






Af-l 



(bq 



2\-l 



{bq 





Af-1\-1 



(4.38) 



(4.39) 



(l^W|f/ = (l,f/oi,--- ,Uo,N-l 



10) 



(4.40) 
(4.41) 
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5 Representation: the case where q = I 



In this section a representation for g = 1 is given. We consider matrices D and E instead 
of d and e. Taking the hmit g — 1 in the representation in section 14.21 yields 



D 



^0 ^0 







D'o D\ D\ 
D\ D\ 



El 







K 


E\ 


E[ 





E\ 


E\ 



(5.1) 



(iy| = (i,o,o,---), |r) = (i,o,o,---f, 



(5.2) 



where 



E>1 



a + 5 + n{aP + 2aS + jS) 



/? + 7 + n{aP + 2/^7 + 7^) 



a 



a + 7 



(« + 7)(/? + 5) 
[in + l){\ + n + l)f\ El 



1 



a + 7 



(« + 7)(/5 + 5) 

\{n + l){\ + n + l)f^ 



[{n + 1)(A + n + 1)]^/' , El = [{n + 1)(A + n + 1)]^/' , (5.3) 



P + 6 



with 



The eigenvector of C is 



a + P + -f + 5 _ 
"(a + 7)(/3 + 5) ■ 



where the elements are 



C(x) = (-1)" 



n!r(A + l) 
_r(A + n + l) 



1/2 



(5.4) 



(5.5) 



(5.6) 



The orthogonal polynomials in this case are the Laguerre polynomials L^n\x). The recur- 
rence relation of lI'^^ (x) is 

{n + l)L^^Ux) - (2n + A + 1 - x)L^^\x) + {n + X)Li%{x) = 0, (5.7) 

with 4^^(x) = 1 and lL^^(x) = 0. In our notation it is rewritten as 

C\i{x)) =x\i{x)). (5.8) 

The orthogonality relation reads 

^dx x^e-^Ll^\x)Ll^\x) = ^^""^^^^^ mn, (5.9) 
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for A > — 1. It is rewritten as 

1 = /i-^ / dx a;V^|£(x))(£(x)|, 



where /z = r(A + 1). We can obtain the form factor as 



where 



{£{x) \D\i{y)) = fi-^ (ci + C2x)x-^e''6{y - x) + C3 hm J{x, y, z) 
L z— +1 



Cl 



a/5 — 7(5 



(a + 7P(/3 + 5)^ 



C2 



a + 7 



C3 



a/? — 75 



(a + 7)(/3 + <5)^ 



(5.10) 



(5.11: 



(5.12) 



1 - z 



[xyzlx+i{x,y,z) - yzlx{x,y, z)], 



Ix{x,y,z) = 



r(n + A + 1 
exp 



-L^^\x)L^\y)z- 



n=0 

{xyz)'^/"^ 



I- z 



Here I\{z) is the modified Bessel function. 



{x + y)z\ f 2{xyzy/'^ 



1 - z 



(5.13) 



(5.14) 



6 Partition Function Zi and Bulk Quantities 

Applying the formulae described in the previous sections, the partition functions of the 
ASEP, Zl, are obtained in the integral form. The results are of course free from the choice 
of representation. 



6.1 The case where q ^ 1 

We consider the representation in section 14.21 since the boundary parameters are unre- 
stricted. First, consider the g < 1 case. Inserting the identity operator ()4.27|) in ()2.4p . 
then we have 

Zl = (WlC'lV) 
dz 



ho' {{z + z-')/2) {W\p {{z + z-')/2)) {p {{z + z-')/2) \C^\V) 

ho' i mP ((^ + ^"')/2)) (P {{z + z-')l2) \V) 



(J Attiz 
dz 



'l + z){l + z- 



1-q 



c 



Aniz 



w {{z + z-^)/2) 



;i + ^)(i + ^- 

1-g 
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In the third hne, we used fj4.26|) . The integral formula of the partition function is now 
obtained: 



(J 4:Triz {az, a/z, bz, b/z, cz, c/z, dz, d/z] q)oo 



(6.1) 



where the integral contour C is a closed path which encloses the poles aX z = aq^, bq^ , 
cq^, dq^{k G Z+) and excludes the poles aX z = {aq'^)~^, {bq^)~^-, {cq^)~^-, {dq^)~^{k G Z+). 
The partition function with the fugacity ()2.12|1 is similarly obtained with replacements 
a a^'^i b bC,, c — > c^~^, d dC,: 



dz 



{z\z~^;qUil + ^z){l + ^z-')/{l-q)Y 



Attiz {a^ ^z,a^ ^ /z,bC,z,bC,/z,cC, ^z,cC, ^/z,dC,z,dC,/z;q)oo 



(6.2) 



The integral is hard to evaluate, but its asymptotic behavior in the thermodynamic 
limit L ^ 1 is obtained without difficulty. We should consider three cases: (A) a > 1 and 
a > b, (B) 6 > 1 and b > a and (C) a < 1 and b < 1. 

(A) a > 1 and a > b: In the integral, the dominant contribution comes from the residue 
of the poles at 2; = a,a^^ of order [(1 + a)(l + a~^)/(l — q)]^. Therefore, we have an 
approximation of as 



(a 2;g)^ 



{q,ab,ac,ad,a ^b,a ^c,a ^d;q)c 



;i + a)(l + a 



Similarly, for the partition function with fugacity, we have 



(a-'e-, q) 



(q,ab,ac^ '^,ad,a ^b^'^,a ^c,a ^d^'^^q), 



l-q 



Then, the bulk quantities ()2.11|) . ()2.13p and ()2.14|) are calculated as 



J^(l-g)- 



(B) b > 1 and b > a: The results are similarly obtained as (A): 



{b-'-q\ 



{q,ba,bc,bd,b ^a,b ^c,b ^d;q)c 

(fe-'r^g)oo 



l-q 



-1^ 



(g, ba, be, bd^'^, b~^a^~'^, b'^c^'"^, b~^d; q)^ 



[i+be){i+ b 

l-q 



Then, we have 

J~ (1 



9) 



{l + bf 



(P) 



1 + 6' 



[l + bfL 



(6.3) 

(6.4) 
(6.5) 

(6.6) 
(6.7) 

(6.8) 



(C) a < 1 and b <1: The contour integral is just over a unit circle and can be approximated 
by the steepest descent method: 



[a,b,c,d;q)l^y/7TL^/^ \1 - q 
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(6.9) 



Then, we have 



2 

oo 



[(i+0(i+r^)] 



3/2 



2 

oo 



J 



1-g 

4 ' 



(P) - ^, 



20FL3/2 



l-q 



.10) 



(6.11) 



The results for g > 1 are directly from the above expressions considering the reflection 

symmetry (j4.4j) . 

(A') > 1 and b'^ > a'^: 



J 



a ^ > 1 and a 



-1 



> b 



(l + &)2' 



-1. 



(P) 



;i + &)2L- 



(6.12) 



J 



(1 + a) 



(P) 



1 + a' 



(Ap^) 



(l + a)2L 



(C 



< 1 and 6 1 < 1: 



J 



1-9 



-1 



(P) - 2' 



(6.13) 



(6.14) 



The minus sign in the front of the current is due to the parity change under the reflection. 

In this way, we have three phases altogether. The phase A is called the low-density 
phase, the phase B the high-density phase, and the phase C the maximal-current phase. 
In fact, the average of particle density increases in the order of phase A, phase C, phase 
B. The particle current takes the maximum value in phase C. Figure El shows the phase 
diagram of J and (p). The same diagram as that in is reproduced. We observe that 
across the line between the phase A and B, the first derivative of J is discontinuous and 
(p) itself is discontinuous. Across the lines between phase A and C, or phase B and C, 
the second derivative of J is discontinuous and the first derivative of (p) is discontinuous. 
Notice that J is given by the partition function Zl and (p) is given by the first derivative 
of Zl- Then, in the terminology of statistical mechanics, we can say the transition is of 
first order across the line between phase A and B, and is of second order across the line 
between phase A and C, or phase B and C. It is interesting to discuss some extreme limit 
cases. When q > 1 and S = 0, or equivalently d = 0, the situation is that there is no 
particle supply at the right boundary. Since the particles are likely to move leftward in 
the bulk and can move out the system at the left boundary, the system becomes short of 
particles in the end. In this extreme remarked in the end of section El the weight 

function in the integral should be changed into the form of (j3.24j) . The partition function 
Zl is then 



du (iaq ^e^^—iaq ^,ibq ^e^,—ibq ^,icq ^e^,—icq g 
logg 



-q ^e^", -q q i)^ 



X 



2 + 2i sinh u 



(6.15) 
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Similarly, when g > 1 and 7 = 0, or equivalently c = 0, the situation is that there is no 
escape of particles at the left boundary. More and more particles accumulate near the 
left "blocking" boundary and finally the system is stopped with jammed particles. The 
partition function in this case is almost the same as the above. Just replace c by d. We 
remark that the case where q> 1 and c = = (or, 7 = 5 = 0) was treated in ^7] . 

6.2 The case where q = l 

The partition function for g = 1 is given by 

1 



r(A + l) 
r(A + L + l 

r(A + i) 



I (6.16) 



The current is calculated as 



Instead of the average density ()2.13p . we compute exactly the particle density profile (rfc) 
in the next section. 



7 The n-Point Function 

The n-point functions are obtained in the integral form as well. In general, we find that 
the n-point function contains n + 1-multiple integrals. 



7.1 The case where q ^ 1 

Inserting the identity operator ()4.27p n + 1 times altogether between D's in the matrix 
products expression and using the form factor ()4.32|) . we obtain the n-point function for 
ii < ■ • • < jn, 



1 



lim 



;i -g)^Ziei, -,en-i 



■71+1 

n 

,m=l 



dZr, 



Cm ^T^i^m 



jm jm — 1 1 



n+1 
m=l 

{azi, a/zi, czi, c/zi, hZn+l-, b/Zn+l, dZn+l, d/zn+i;q) c 

1 



n 

m=l 



"'^(^m+1 ^m) 



(^m' ^rr^J Q)cx> 

cx3(^m+l ~l~ 1/^m+l ^mZ-m ^m/ ^m) 



(7.1) 
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where jo = 0, jn+i = L + 1 in the product. Regularization parameters ^m's are included 
since otherwise this integral becomes a singular integral. The contours Cm(m = 2, ■ ■ ■ ,n + 
1) encloses the poles at Zm = ^m-iZm-i<l'' , im-iz^-iQ^ (k E Z+) and excludes the poles at 
Zm = {^m-iZm-i<l'')~^ , ('Cm-i-Zm-i'?'^)^"'^ ^ ^+)- addition, for m = 1, Ci encloses the 
poles at Zi = aq^, cq^ {k G Z+) and excludes the poles at Zi = {aq^)~^, {cq'')~^ {k G Z+). 
Also, Cn+i encloses the poles at Zn+i = hq^, dq^ {k G Z+) and excludes the poles at 
Zn+i = {bq^)~^, {dq^)"^ {k G Z+). Note that the result should be independent of the order 
of inserting of the identity operator. The above choice of integral contours is when the 
insertions are done from left to right. 

The results for the g > 1 case are simply obtained from the reflection symmetry ()4.4|) . 

7.2 The case where q= 1 

Applying formulae in sectional we also obtain the n-point function for q = 1: 



1 



n+l 

,A „-Xr. 



m=l 



n 



72+1 

:jm—jm — l — 1 



m=l 



U)] , (7.2) 



X 

m=l 



where jo = 0, = L + 1 in the product. In particular, for the one-point function the 
integral can be performed directly. Thus, we obtain 

(rO ^ ^ - "^-J' ^ (-^ + k-l). (7.3) 



a + 7 X + L{a + -f){f3 + 6) \a + 
for k = 1, - ■ ■ ,L. This shows a linear profile of particle density. 



8 Conclusion 

We have studied the steady state of the ASEP for the most general case of open boundary 
conditions where particles hop in both directions and are injected and ejected at both 
boundaries. Finding out an explicit representation of the algebraic relations related to 
the Askey- Wilson polynomials, we have obtained the partition function in the form of the 
moments with respect to the weight function of the Askey- Wilson polynomials. This is 
an extension of the previous work for two degrees of freedom to the most general case. 
Furthermore, we also have found a novel integral formula for the ra-point functions. 

From the asymptotic behaviors of the partition function, the current in the thermody- 
namic limit has been calculated. The phase diagram is given and the boundary-induced 
phase transitions are observed. For the n-point functions, a limit procedure is included in 
the multiple integrals. The asymptotic evaluation of them is more difficult, but could be 
done by generalizing the methods in [TH] . 

The advantage of our method using the orthogonal polynomials is that, once the corre- 
sponding polynomials are known, one can use the whole storage of knowledge on them. It 
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should be stressed that the Askey- Wilson polynomials are considered to be one of the most 
important orthogonal polynomials because all important classical orthogonal polynomials 
are obtained as special/limiting cases of the Askey- Wilson polynomials. 

The discovery of the connection to the Askey- Wilson polynomials is yet another evi- 
dence of a deep mathematical structure of the ASEP, which has been already suggested in 
several previous works. For instance, the ASEP can be considered as an exactly solvable 
spin system and the Bethe ansatz methods has been applied successfully [30, 31J. The 
connection of the matrix method to the Zamolodchikov-Faddeev algebra was pointed out 
in [S2]- In addition, more recently, fluctuations of the current of the ASEP have been 
studied using the connection to the random matrix theory jHSl ■ would be of great 
interest to further clarify the interrelationship among these fascinating connections. 
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Figure 1: The ASEP with open boundaries 
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Figure 2: The phase diagram of the ASEP. There are three phases; phase A: low- 
density phase, phase B: high-density phase and phase C: maximal-current phase, (x, y) = 
(a~^,6~^) if g < 1, and {x,y) = {b,a) if g > 1. a and b are defined in (jHH). The double 
line is for the transition of the first order and the solid line is for that of the second order. 
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